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Abstract 

Quantum critical behavior in 2+1 dimensions is established via holographic meth- 
ods in a 5+1-dimensional Einstein gravity theory with gauge potential form fields of 
rank 1 and 2. These fields are coupled to one another via a tri-linear Chern-Simons 
term with strength k. The quantum phase transition is physically driven by the ex- 
pulsion of the electric charge from inside the black brane horizon to the outside, where 
it gets carried by the gauge fields which acquire charge thanks to the Chern-Simons 
interaction. At a critical value k = k c , zero temperature, and any finite value of the 
magnetic field, the IR behavior is governed by a near-horizon Lifshitz geometry. The 
associated dynamical scaling exponent depends on the magnetic field. For k < k c , the 
flow towards low temperature is governed by a Reissner-Nordstrom-like black brane 
whose charge and entropy density are non-vanishing at zero temperature. For k > k c , 
the IR flow is towards the purely magnetic brane in AdS6- Its near-horizon geometry 
is AdS4 x R 2 , so that the entropy density vanishes quadratically with temperature, 
and all charge is carried by the gauge fields outside of the horizon. 



1 This work was supported in part by NSF grant PHY-07-57702. 



1 Introduction 



Via the AdS/CFT correspondence, classical gravity can be used to model strongly interact- 
ing matter at finite temperature and finite charge density. Much activity in recent years has 
been devoted to finding classical solutions to various theories of gravity whose behavior bears 
some relation to interesting condensed matter systems. To the extent that strong coupling 
dynamics plays an important role in a given condensed matter system, the AdS/CFT cor- 
respondence is one of the few computational tools available. Holographic computations of 
thermodynamics and transport are relatively easy to carry out, since they reduce to questions 
involving classical gravity. See [U |2] for reviews. 

Quantum critical systems are particularly interesting to study in this regard. On the 
one hand, quantum criticality appears to be an important feature of various real world 
systems [3]. On the other hand, the inherent IR universality of such quantum critical points 
implies that it is meaningful to model these systems by the kinds of supersymmetric gauge 
theories that naturally arise in the AdS/CFT correspondence, even though their UV details 
may bear no relation to any realistic condensed matter system. 

Experimentally, quantum critical points are found by tuning various control parameters, 
such as a pressure or a magnetic field, and one can attempt to map out a phase diagram 
in terms of these parameters. On the gravity side, it is possible to proceed similarly. One 
important driver of holographic quantum criticality is related to how charge is carried by 
such systems; in particular, the qualitative nature of the charge carriers can change under a 
variation of the control parameters, and this can signal a quantum critical point. 

In the AdS/CFT correspondence, the density of charge in the boundary theory is equal 
to the amount of bulk electric flux that pierces the boundary. In the simplest bulk solution, 
the electrically charged Reissner-Nordstrom (RN) black brane, Gauss' law implies that these 
flux lines must all emanate from the horizon of the black brane, and so the charge carriers 
are hidden behind the horizon. The physical interpretation in the boundary theory is poorly 
understood, partly due to the fact that this RN solution has a large ground state entropy 
density]^] An alternative scenario is for there to be a source on the right hand side of Gauss' 
law, in which case the electric flux lines can be generated in the bulk as one moves out 
radially. One way this can occur is if there is a condensate of a charged scalar field in the 
bulk; this spontaneously breaks the U(l) gauge symmetry and corresponds on the boundary 
to a holographic superfluid phase [5j. A second possibility is for there to exist a sea of 
charged particles in the bulk; solutions of this type have been studied in [6j. 



In the present paper, we are primarily interested in yet another possibility, in which 
Gauss' law is modified by Chern-Simons interactions. In the presence of a Chern-Simons 

2 It has been suggested, e.g. in [3], that one should think of these as fractionalized phases of matter, in 
analogy with the fact that horizons are known to describe the deconfined phases of gauge theories in the 
AdS/CFT correspondence. 
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term, gauge fluxes can themselves carry charge, and this allows for charge to be located 
outside an event horizon in the absence of explicit charged particles, and without breaking 
the U(l) symmetry. This mechanism played a key role in the magnetic quantum critical 
points studied in El |H| (see [TO] for another example of magnetic quantum criticality in a 
holographic setting). There, the bulk theory was five-dimensional Einstein-Maxwell-Chern- 
Simons theory. Solutions were obtained that are dual to a boundary gauge theory at finite 
charge density and external magnetic field. In the bulk, the simultaneous presence of electric 
and magnetic fluxes activated the Chern-Simons term, allowing bulk charge to be carried in 
the manner described above. The strength of this effect was varied by changing either the 
Chern-Simons level k or the value of the magnetic field B, and a rich phase structure was 
uncovered. Most interestingly, for sufficiently large k, a magnetically tuned quantum critical 
point with nontrivial scaling properties was discovered, resembling real world systems in a 
number of ways [7j. 

The critical theories studied in [TJE] are effectively 1+1 dimensional: while the underlying 
gauge theory is defined in D=3+l, at the critical point massless propagation takes place 
along a single spatial direction. In the bulk, this was seen from the existence of a near- 
horizon geometry of the form WAdSs x R 2 , where WAdS3 is a "null warped" deformation of 
AdS3 space [HJ IT2] . The R 2 was supported by the nonzero magnetic field. While quantum 
criticality in 1+1 dimensions is certainly important, it is also interesting to generalize to 
higher dimensions. Particular interest attaches to the case of 2+1 dimensions, which in the 
real world is relevant to various layered materials such as the cuprates. Our goal in the 
present work is thus to find a magnetic quantum critical point in 2+1 dimensions, driven by 
the presence of Chern-Simons terms. 

1.1 Overview of the set-up 

The simplest setup exhibiting the desired behavior appears to be six-dimensional gravity 
coupled to a one-form potential A and a two-form potential C. Assuming gauge invariance 
under A — > A + dA^ and C — > C + dAc, the most general two-derivative action is, 

S = J d 6 x^ (r~^ + F mn F mn + ±G MNP G MNP ^j + S cs + S hndy (1.1) 

with the Chern-Simons (CS) term Scs given by, 

s - = i4/ CAFAF (L2) 

Here F = dA and G = dC are the two-form and three-form field strengths respectively. 
While such an action can plausibly arise within string theory, for the time being we will 
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regard this as the result of a bottom-up construction]^] Although the dual boundary theory 
is defined in 4+1 dimensions]^] as discussed above our primary interest is in finding solutions 
corresponding to quantum criticality in 2+1 dimensions. This is achieved by introducing a 
magnetic field, which acts to freeze out two of the spatial dimensions in the IR. 

To turn on a nonzero charge density p and magnetic field B we consider a two-form field 
strength of the form 

F = E(r) dr A dt + Bdx 3 A dx 4 (1.3) 

The value of the charge density p is read off from the asymptotic behavior of E, namely 
p ~ r 4 E(r) as r — > oo. The CS term induces a nonzero three-form field strength, 

G = G 1 (r)dr A dx 1 A dx 2 (1.4) 

The full solution is taken to be translationally invariant along all boundary directions, as 
well as rotationally invariant in the x 1,2 and x 3 ' 4 planes. The form- field equations of motion 
exhibit a source in Gauss' law due to the CS term, 

d* F — 2kF A G = 
d* G + kF A F = (1.5) 

In particular, the product BG\ acts as a charge density sourcing E. 



1.2 Summary of the results 

Given this setup, we can try to characterize the zero temperature solutions of this theory as 
a function of the charge density p, magnetic field B, and Chern-Simons level k. Actually, due 
to scale invariance there are really only two parameters to vary: the dimensionless parameter 
k and the dimensionless ratio, 

B= B - (1.6) 

One way to summarize our results is to give the low temperature behavior of the entropy 
density s for various values of B and k. Let us first note two special cases that are easy to 
understand. For B = the solution is just the electric RN black brane, with a finite entropy 
density at zero temperature. The other extreme is B = oo, attained with E = p = 0. In 



3 A related theory with two-form and three-form potentials and Chern-Simons terms is obtained from 
massive IIA supergravity [15] , and admits solutions with a four-dimensional Lifshitz factor |16l 117] , similar 
to what we find below. 

4 CFTs in 4+1 dimensions are rare, but not unheard of |18j . 
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this case we find a six-dimensional version of the magnetic brane solution studied in [13] for 
five dimensions. The solution develops a near-horizon AdS4 x R 2 geometry which connects 
to the asymptotic AdSg. In both of these cases the three-form field strength G vanishes 
identically, and the Chern-Simons term plays no role. 

For finite nonzero values of B we find that the qualitative behavior depends on the 
magnitude of k relative to its critical value k c — 1/ y/3. We discuss the three cases in turn. 

1.2.1 The case k < k c 

Solutions in this regime can be thought of as deformed versions of the purely electric RN 
solution. At B = 0, we have precisely the RN solution, as G\ vanishes and the CS term 
plays no role. At zero temperature, the near-horizon region becomes AdS 2 x R 4 , and the 
zero temperature entropy density is proportional to the charge density: s oc p. For this RN 
solution, all of the charge is hidden behind the horizon. 

Turning on nonzero B activates G\ and the CS term. At zero temperature the near- 
horizon region again exhibits an AdS 2 x R 4 , but now only a finite fraction of the charge 
is hidden behind the horizon, with the remainder carried outside the horizon by the flux 
term BG\. As k and/or B is increased, this effect becomes more pronounced, with more 
charge transferred from behind the horizon to outside, as illustrated in figure [Tj The transfer 
becomes complete as k reaches its critical value k c , at which point there is a quantum phase 
transition. For k < k c , all of these solutions have a nonzero ground state entropy density. 




Figure 1: As k and B are increased (indicated by the blue arrow), the charge inside the 
horizon of the black hole gets expelled and carried by the gauge fields outside the horizon. 
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1.2.2 The case k > k c 

Here the solutions are generalizations of the purely magnetic brane solution obtained for 
E = G\ = 0. The purely magnetic solution interpolates between a near-horizon A0IS4 x R 2 
region and AdS6 asymptotically, and is the analog of the 5D magnetic brane solution studied 
in [13]. The low temperature entropy density behaves like that of a D=2+l Lorentz invariant 
CFT, and exhibits the expected s oc T 2 . At small but nonzero temperature, the near-horizon 
region has AdS4 replaced by an AdS4 Schwarzschild black brane. 

Turning on nonzero charge density changes the details, but the qualitative picture re- 
mains. The zero temperature solutions always exhibit a near-horizon AdS4 x R 2 region and 
the entropy density is quadratic in the temperature. The dependence on B and k can be 
seen in the behavior of the prefactor, 

s~A{k,B)f 2 (1.7) 

where we now work in te rms of dimensionless quantities, referred to with hats, and to be 



defined precisely in (3.29). It is interesting to explore the behavior of A(k,B) as B tends 
to 0, since precisely at B = (obtained for B = and p ^ 0) we arrive back at the purely 
electric RN solution with its finite entropy density ground state. The coefficient A(k, B) 
must therefore exhibit a singularity, and indeed we find the following behavior, as B — > 0, 

A(A;,5)~c(A;)exp{^} (1.8) 

for real positive functions c(k) and d(k). We present clear numerical evidence for the simple 
but distinctive analytical behavior of A(k, B) (and for d(k) as explained below), even though 
at present the corresponding analytical derivations are not available. 

For k > k c we therefore find that the ground state entropy density of the purely electric 
RN solution is removed by the presence of any nonzero magnetic field B. For any finite B, the 
entropy density will always go to zero quadratically in T as T — )• 0. The finite ground state 
entropy density of the RN solution thus requires a fine tuning, since it is removed by turning 
on an arbitrarily small magnetic field. It is instructive to compare this to the behavior of the 
analogous charged magnetic brane solutions in AdS 5 studied in [H]. There, we conjectured, 
based on preliminary numerical analysis, that the AdSs RN ground state entropy density 
would similarly be removed by the presence of any finite magnetic field. However, subsequent 
high-precision numerical work [7] and analytical calculations [8] revealed that the magnetic 
field needed to be larger than a value B c to achieve this. So in this sense the effect of the 
magnetic field appears to be stronger in the AdS@ compared to AdS 5 . 

The function A(k,B) is also singular as k approaches k c from above. Numerical analysis 
shows that the coefficient d(k) behaves as 

<*(*) - jJ^2 (1-9) 



6 



to remarkably good accuracy in the range k c < k < 1.2. Here, do is a constant of order 1. 
This again signals a quantum phase transition at k c , to which we now turn. 



1.2.3 The case k = k c 

Numerical analysis at k = k c indicates that the low temperature entropy density scales to 
zero as a S-dependent power of temperature, 

{oo as B — > 
(1.10) 
1 as B — > oo 

A plot of z(B) will be shown in figure]^] Such a behavior is characteristic of a 2+1 dimensional 
scale invariant theory with dynamical critical exponent z = z(B). Indeed, precisely at k = k c 
the field equations admit a new near-horizon geometry of the form Lif 4 x R 2 , where Lif 4 is 
a 3+1 dimensional Lifshitz solution [Tj5] with dynamical exponent z. Explicitly, 



ds 2 = —-uo^de + r^Hidx^ + idx^ + Udx^f + idxy) 
F = q dr A dt + b dx 3 A dx 

G = --^r 2tz - x dr Adx 1 Adx 2 (1.11) 

KOqZ 

In terms of z restricted to the range 1 < z < oo, as is required for all fields to be real valued, 
the values of go and bo are fixed as 

_ 10z(z - 1) _ 30, 

The parameters go and bo represent the near-horizon value of the electric charge and magnetic 
field respectively, and are related to asymptotic quantities measured at the AdS6 boundary 
by the full interpolating solution. There is thus a nontrivial expression for B in terms of z, 
or equivalently z(B), whose form is only known numerically; see figure ^ below. 

Numerical analysis shows that as the temperature is lowered at k = k c a near-horizon 



region of the form (1.11) indeed develops and controls the low temperature thermodynamics. 
Also, as B varies from to oo, the dynamical exponent z smoothly varies between oo and 
1, so that the family of Lifshitz geometries interpolates between AdS2 X R 4 and AdS 4 x R 2 . 

In figure[2j we display representative numerical data that illustrate the points made above. 
In the left panel (b, q) denote the values of the electric and magnetic fields at the horizon 
(in a particular coordinate system described in the text). By varying the temperature for a 
fixed value of B we obtain curves in (b, q) space. The curves emanate from the origin, which 
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Figure 2: Left panel: Entropy density s versus temperature T for selected values of k. Right 
panel: Flows towards low temperature of the horizon data q and b, for fixed B = 1, for 
k = 0.1 (red dots), k = k c (black dots), and k = 0.7 (magenta dots). 



corresponds to the high temperature regime. As the temperature is lowered, the curves either 
flow to the electric fixed point at b = for k < k c ; to the magnetic fixed point at q = for 
k > k c ] or to a point on a critical curve with nonzero q and b for k = k c . In the right panel 
we show the corresponding behavior of the entropy density as a function of temperature. As 
T — > 0, the entropy density either goes to a finite value (k < k c ); to zero quadratically in 
temperature (k > k c ); or to zero with a B dependent power between and 2 (k = k c ). 

Examining the overall picture, summarized in figure [3| the story is quite similar to that 
in the AdSs case, except that the roles of k and B are, in a sense, interchanged. In the AdSs 
case a critical point was reached by tuning B to B c , and at the critical point the dynamical 
exponent z depended on k. But in the AdS6 case studied here, the critical point is reached 
by tuning k, and at the critical point the dynamical exponent varies with B. The other 
major distinction is, of course, that for AdSs the critical theory was 1+1 dimensional, while 
in the AdSg case it is 2+1 dimensional. 

A model with related behavior was studied in (20]. There also, transitions were found be- 
tween finite and zero entropy density states at zero temperature, mediated by the appearance 
of a Lifshitz near horizon geometry. 



8 




Figure 3: Phase diagram of entropy density s versus k for all finite B. 



1.3 Organization 



The remainder of this paper is organized as follows. In section [2j we present a general discus- 
sion of candidate theories for holographic quantum critical behavior in 2+1 dimensions, built 
on Einstein gravity plus gauge fields of various ranks. The theory of gauge field potentials 
of rank 1 and 2 in six space-time dimensions with a Chern-Simons coupling, considered in 
this paper, emerges as one of the simplest such models. In section |3j a detailed discussion 
of the action, field equations, symmetries, suitable Ansatz, and reduced field equations for 
the study of thermodynamic quantities is presented. In section |4| the special solutions cor- 
responding to the well-known AdS6 Reissner-Nordstrom black brane, and the novel purely 
magnetic brane in AdS6 are derived. The critical value of k = k c is inferred from the 
behavior of small fluctuations around the purely magnetic brane. The existence of a Lif- 
shitz near-horizon geometry at k = k c is shown in section |5j its extension to a full fledged 
six-dimensional solution is derived numerically, while an analytical derivation of the small 
fluctuation spectrum around Lifshitz is deferred to Appendix [A} In section [6j a detailed dis- 
cussion is presented of our numerical results, and their matching with analytical calculations 
when available (with further details explained in Appendix B). Finally, the quantum critical 
behavior of the system is discussed in section [7j 
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2 Candidates with -p-form fields 



The holographic model for quantum critical behavior in 1+1 dimensions based on 4+1- 
dimensional Einstein-Maxwell-Chern-Simons [H [8] draws on a number of fundamental dy- 
namical properties. 

1. The presence of a constant background magnetic field induces an IR flow towards a 
lower-dimensional AdS theory; 

2. The gauge field can carry charge thanks to the Chern-Simons interaction. This mech- 
anism permits charge to be expelled from the inside of the black brane horizon to the 
outside where it is carried by the gauge field. 

3. Tri-linearity of the Chern-Simons action supplies the IR theory with an effective Chern- 
Simons theory for the low-energy fields. 

The goal of this paper is to generalize and extend this construction to space-times where 
quantum criticality takes place in field theories of dimension higher than 1+1. To do so, we 
shall want to respect the basic tenets above. Given these assumptions, which theories can 
we construct, while maintaining the basic assumption that no scalar fields be present? 

2.1 Gravity models with a single form- field 

The Einstein-Maxwell-Chern-Simons theory of [jj |8], formulated in terms of a single 1-form 
field potential, admits a straightforward generalization to a theory with a single p-form field 
C whose field strength form is denoted by G = dC . The corresponding action is given by, 

S = -tttV / d ° x ^9 (R — A+ \G\ 2 ) + S cs (2.1) 
where A is the cosmological constant, Scs the Chern-Simons term given by, 

s - = dfe/ CAGAG < 2 - 2 > 

and the boundary action has been suppressed. Matching of dimensions requires D = 3p + 2. 
The Einstein-Maxwell-Chern-Simons theory of 0, |8] corresponds to p — 1. When p > 1, we 
assume that G has a boundary magnetic field that generates a flow from AdS£>=AdS3 P+ 2 in 
the UV to AdSi)_p_i=AdS2p+i in the IR. For p = 2, we have D = 8 and a flow towards AdSs, 
while for p = 3, we have D = 11 and flow towards AdS7. The latter case might be of special 
interest as arising from M-theory on a S 4 sphere with magnetic field. All of these schemes 
may be interesting, but still they lack the presence of the fundamental electro-magnetic field. 
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2.2 Gravity models with two form- fields 



Gravity models with two different form potential fields A and C, respectively of ranks p' 
and p, and field strengths F = dA and G = dC, offer a further extension of the Einstein- 
Maxwell-Chern-Simons models of [7, 8\r \ The general action is as follows, 



(2.3) 



S = -^~q^ J d ° x V9 (R ~ A + \F? + \G\ 2 ) + S cs 
where A is the cosmo logical constant. All possible tri-linear Chern-Simons terms are, 



AAA 
CS 



qAAC 
°CS 



qACC 
°CS 



<1 



ccc 

CS 



AttG d 
AnG D 

h 

4irG D 
AnG D 



A A F A F 
A A F AG 
A AG AG 
C AG AG 



D = 3p + 2 
D = 2p' + p + 2 
D = p' + 2p + 2 
D = 3p + 2 



(2.4) 



In the presence of two different form fields of the same rank p = p', all of the above Chern- 
Simons terms are allowed for D = 3p + 2, and the theory is characterized by 4 independent 
dimensionless couplings hi, k 2 , k 3 , k<±. On the other hand, if p ^ p', and we take p' < p, then 
for given D only a single one of the above Chern-Simons couplings can occur. 

Assuming that one of the potentials A is electro- magnetic, we set p ! = 1, and p > 1. For 
p = 1 we recover the D = 5 theory considered in [3 E] except now with two gauge fields and 
a mixed Chern-Simons term; theories of this type were studied in [221 123]- For D > 5, each 
value of p > 2 yields a unique theory with Chern-Simons term S^ c in dimension D = 4+p, 
and a unique theory with Chern-Simons term S^ c with dimension D = 3 + 2p. The next 
simplest case has p = 2, D 
here. Another interesting case is p = 2, D 



6 and Chern-Simons term S^g and is the one considered 



7 and Chern-Simons term S^g ° . 



3 The Einstein-Maxwell-two-form-field system 

In this section, we present the Einstein-Maxwell-two-form-field system in detail, including 
the action and field equations. To study thermodynamics in the presence of a uniform 
background magnetic field and electric charge density, we focus on solutions which are in- 
variant under translations, and certain rotations, of the boundary theory. We obtain the 
corresponding Ansatze and reduced field equations. 

We will be assuming invariance (up to boundary terms) of the action under the gauge invariance A — > 
A + cIAa, C — > C + dAc, although further generalizations could be envisaged as well. 



ff 



3.1 Field contents and field equations 



The fundamental fields of the 6-dimensional theory are: the metric QmNi the Maxwell field 
Am, and the 2-form field Cmn = —Cnm- The action consists of the Einstein-Hilbert term, 
standard kinetic terms for the fields A and C, and the only possible Chern-Simons term 
afforded by this structured 



S = - 167 ^r / d e x^g (r~§ + F mn F mn + \g mnp G mnp ^ + S cs + S hndy (3.1) 
with the Chern-Simons term Scs given by, 

The explicit form of the boundary action SWay will not be needed. Here, Gq is Newton's 
constant in 6 dimensions, L sets the scale for the cosmological constant, and k is the sole 



remaining dimensionless coupling. The action (3.1) is invariant under the gauge invariance 
A — > A + dAji, C — > C + dAc, up to a boundary term generated by Scs] in fact, (3.1) is 
the most general two-derivative action with this property. The field strength and potential 
forms, and their components, are related as follows, 

A = A M dx M F = dA = ^F MN dx M A dx N 

C = \c M Ndx M A dx N G = dC = -Gmnp dx M A dx N A dx p (3.3) 

2 6 

Einstein's equations may be cast in the following form, 

Rmn = —2F MP F N P — G M pqG n pq + g M N (j^ + ^F P qF pq + -G pQ R G PQpS j (3.4) 

Note that the contribution of the F-field is traceless in 4 dimensions while that of the G-field 
is traceless in 6 dimensions, as we should expect. The Maxwell equations in the presence of 
the 2-form field C are given by the Bianchi identity dF = 0, together with the field equation, 

VmF mn = e MNP( * RS F PQ d M C RS (3.5) 

^ V 9 



6 Throughout, Einstein indices range over M,N,P — 0,1,2,3,4,5, and we define g = —det(gnjN)- Iso- 
lating the holographic radial direction r, and the time direction 0, we adopt the following convention for 
the orientation, e r oi234 = 1- Expressed in Einstein indices, the dual *G of a p-form G with components 
G Ml -M p , is defined by {*G) Nl ... Ne _ p = £m 1 -m p n 1 ---n^ p G Mi '" m " . 
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The 2-form field equations in the presence of the Maxwell field are given by the Bianchi 
identity dG = 0, together with the field equation, 

V P G PMN = -±-e MNPQRS F PQ F RS (3.6) 
4 v 9 

In form notation, these equations become, 

d* F — 2kF A G = 

d* G + kF A F = (3.7) 

Henceforth, we set L — 1, and refer to the fields F, G as the gauge fields of the system, and 



to (3.7) as the gauge Geld equations. 



3.2 Ansatz with translation and rotation symmetries 

The precise symmetry requirements will depend upon the physical questions addressed. 
Throughout, we shall be interested in thermodynamic quantities (postponing problems in- 
volving correlators to later work) evaluated in the presence of a uniform magnetic -F 34 = B 
field in the x 3 ,x 4 directions, as well as a physical electric charge density p. Therefore, the 
natural symmetries to be imposed will always include the followingj^] 

1. Translations in t, Ob . Ob ■ Ob j X j 

2. Rotations in the x 1 x 2 -plane; 

3. Rotations in the x 3 x 4 -plane. 

The Ansatz for the general metric invariant under the above symmetries is given as follows, 

ds 2 = d ^- Udt 2 + e m ((dx 1 ) 2 + (dx 2 ) 2 ) + e 2V > ((cfe 3 ) 2 + (dx 4 ) 2 ) (3.8) 

This metric is expressed in a specific gauge choice for the holographic coordinate r in which 
9rr9tt — — 1- The functions U, V±, V% depend only on r by translation invariance. Subject to 
the same symmetries, the gauge field Ansatze are given by, 

F = Edr A dt + B x dx x A dx 2 + B 2 dx 3 A dx 4 

G = (pxdr + G 2 dt) A dx 1 A dx 2 + (G 3 dr + G 4 dt) A dx 3 A dx 4 (3.9) 



7 The field equations are also invariant under the sign-reversal of F. We will not be led to impose this 
symmetry on our holographic solutions, since we know that this discrete symmetry will be broken by a finite 
electric charge density. 
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Translation symmetry again restricts all coefficients to be independent of 
To satisfy the Bianchi identities dF = dG = the coefficients Bi, B 2 , G 2 , G4 must be inde- 
pendent of r. The remaining coefficients E, G%, G3 are generally functions of r. 

Throughout, one of the magnetic fields will be kept non-zero, and used to generate the 
RG flow from AdSg to a 2+1-dimensional IR boundary theory The other magnetic field, 
which lives in this 2+1-dimensional IR theory, will be set to zero for simplicity. We choose, 

£i = B 2 = B (3.10) 

where B 7^ 0. (Note that for B\ = B 2 , exact solutions may be found which are closely related 
to the dyonic AdS4 Reissner- Nordstrom black brane |13j.) 

3.3 General reduced gauge field equations 

It is straightforward to compute the reduced field equations for the gauge fields, and we find, 

Ml = (Ee 2Vl+2V2 )' + 2kBG x 

M2 = [dUe 2 ^- 2 ^)' + 2kBE 

M3 = (GzUe 2Vl ~ 2V2 )' (3.11) 

as well as the relation, 

2kBG 2 = (3.12) 

Note that the Bianchi identities were solved already in the preceding subsection. Since we 
are assuming throughout that kB 7^ the field equations imply G 2 = 0. Furthermore, 
equation M3 may be integrated once, to obtain, 

G 3 (r) = _^_ e -^«+2V a (r) (3 13) 

where C3 is an arbitrary integration constant. 

3.4 General reduced Einstein equations 

The reduced Einstein equations were computed using Maple. The equation for the R rt 
component is simply given by, 

2G 1 G 2 e- Wl + G 3 G A e~ m = (3.14) 
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Since we have already deduced the relation G2 = from the reduced gauge field equations, 
we conclude that we must have G3G4 = 0. With the help of (3.11), this condition becomes 
a relation between constants, 



C3 G4 







(3.15) 



Eliminating Gs(r) using equation (3.11), and replacing the Ru equation by the combination 
(R rr — R tt )/(AU), we find the following reduced Einstein equations, 



El 



El 



E3 



E4 



= V{' + (VI) 2 + + {Vtf + G(e 



-4Vi 



+u- 



Cn e 



-4Vi 



+ Gle~^ 



= uv" + u'v; + 2uv[ 2 + 2uv;Vz - 5 



+f/G 2 e -4^ 1 _ [/ -1^2 e - 



l E 2 _ 1^2 -4^ 

2 2 



4Vi 



G\e~^ 



= UV£ + U% + 2UV? + 2UV! 1 V[ - 5 + ^E 2 + ^B 2 e~ 4V2 



- UG l e -^ + U -i ( c 2 e 



-4Vi 



G 2 e" 4V2 



(3.16) 



= U" + 2[/'(V 1 / + Vf) - 10 - 3£ 2 - S 2 e" 4Va 
-2UG\e-^ - 2U- 1 (4 e~ m + G\ e~ 4V *) 

The constraint equation is obtained as the combination E2 + E3 — UE1, and is given by 

con = u'{v{ + v£) + u(y; 2 + k,' 2 + w{v£) - 10 + e 2 + 5 2 e - 4y2 

_^G?e- 4yi - f/" 1 (c 2 e" 4 ^ + G 2 e" 4 ^) (3.17) 

Using Maple, one can check that the derivative of the constraint equation CON vanishes in 
view of the other Einstein and gauge field equations. 



3.5 Reduced equations for solutions with a smooth horizon 



We look for a smooth finite temperature horizon, at which U has a single zero. Other 
solutions will certainly exist, but it seems plausible that only solutions with smooth horizons 



are of physical interest to us. From Einstein equation El in (3.16) it is clear that we need 
C3 = G4 = if we are to assume that V\, V%, and G\ are smooth functions at the horizon. This 



condition, of course, automatically solves the Einstein equation (3.15). With this restriction, 
the gauge fields simplify as follows, 



F = E(r)dr A dt + Bdx 3 A dx 4 
G = Gi(r)dr A dx 1 A dx 2 



(3.1* 
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Equation M3 in (3.11 ) is satisfied automatically, while equations Ml and M2 are unchanged; 



we repeat them here for later convenience. 

Ml = (Ee 2Vl+2V2 )' + 2kBd 

M2 = (GiUe 2V2 - 2Vl )' + 2kBE (3.19) 



With c 3 = G 4 = 0, the reduced Einstein equations of (3.16) simplify and become, 
El = V{' + (V[) 2 + + (V 2 ') 2 + Gie~ m 

E2 = UV" + U'V{ + 2UV! 2 + 2UV[VL - 5 + 1 £ 2 - -SV 4 ^ + C/G?e" 4Vl 

£3 = UV% + C/V 2 ' + 2UVf + 2UV! 1 V{ - 5 + ^E 2 + ^B 2 e' 4V2 - UG\e~ Wl 
EA = U" + 2U'(V{ + Vz) - 10 - 3E 2 - 5 V 4F2 - 2[/G 2 e" 4V ' 1 (3.20) 

while the constraint is given by, 

CON = U'{V; + V£) + U(V[ 2 + + 4^'^) - 10 + E 2 + S 2 e" 4V2 - UG\e~ Wl (3.21) 
3.6 Horizon Data 

For the purpose of carrying out numerical analysis, we need to fix the remaining freedom 
in the choice of coordinates. We put the horizon at r = 0. Using the freedom to rescale 
time t (while compensating by a rescaling of r), as well as x 1 * 2 and x 3 ' 4 , we shall normalize 
the components of the metric by setting, 

U(0) = V 1 (p) = V 2 (p) = Q U\0) = 1 (3.22) 

and use the following notation for the components of the gauge fields, 

B = b E(0) = q G 1 (0)=g (3.23) 

Solving E2, E3, CON, and M2 at the horizon produces expressions for the derivatives 
V{(0), V^O), as well as a relation giving g in terms of the other data, 

g = -2kbq Vl(0) = 5- l -q 2 + h? 

YM = 5 - \q 2 - \b 2 (3.24) 

Thus, the free parameters, as specified at the horizon, are now b and q. We should find 
a unique solution for each choice of (b,q), though the requirement of global regularity will 
restrict the allowed range of these parameters. 
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3.7 Asymptotic data and physical quantities 

For every regular solution, with specified values of (6, q), the asymptotic behavior as r — > oo 
will approach the geometry of the boundary of AdS@. The r — > oo asymptotic behavior of 
the components of the metric takes the familiar AdS6 form, 

Tj(r) e 2Vl ^ p 2 ^W 

~^^v 2 (3.25) 

The asymptotics of the gauge field strengths fix the physical charges of the boundary theory, 

E(r)r 4 ^ qoo Gi(r)r 2 ^ fc (3.26) 

Physical quantities, namely the electric charge density p, the temperature T, the entropy 
density s, the magnetic field B p and the 3-form charge density G p are obtained as follows, 

P = <?oc T = ]- B p = — 

47T V 2 

' G p = 9 — (3.27) 



^V\V 2 Vi 

In our definition of the entropy density we are absorbing Newton's constant: s = (G$S)/Yo\. 
We also comment on the form of p. The electric charge of the system is given by / *F, 
where the integral is over a constant time slice at the boundary. Since F = Edr A dt we 
have -kF ~ V\V 2 r 4 Edx x A dx 2 A dx 3 A dx 4 . Dividing the integrand by ViV 2 , to obtain the 
proper charge density, gives q^. Note that the temperature in these coordinates is fixed at 
the constant value T = l/(47r). Since T is dimensionful we can always scale coordinates 
to bring any nonzero temperature to this value. Of course, this still allows the physically 
relevant dimensionless temperature (denoted below as T) to vary. 

As always, the only physically observable quantities in the conformal boundary theory 
are dimensionless quantities. The above physical quantities all have non-trivial scaling di- 
mensions under r — > £r, and x^ — > x^/£, which are given by, 

p -> tp T £T B p ^ £ 2 B P 

s ->• ts G p ->• £ 3 G P (3.28) 

Thus, only suitable dimensionless ratios can be physical observables. When B p ^ 0, as will 
be almost always the case here, we may use B b to obtain dimensionless physical observables, 

T ^U~ 2 ^ B p b 



T = —= = -*-4= B = 



An^/b yfp v 2 ^fq^ 

v 2 „ _ G p _ g< 



B 2 P ^b 2 y fm Vl V63 



(3.29) 
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Given that the regular supergravity solutions are completely specified by the two parameters 
(b, q) of the near- horizon data (whose domain is governed by regularity restrictions), together 
with the unique dimensionless coupling k, it follows that amongst the physical data T, B, s, 
and g, there must exist two dynamical relations, or equations of state. Choosing T and B 
as independent parameters (as well as k), the dynamics of the system then determines, 

s = s(k,f,B) 

g = g(k,f,B) (3.30) 

These functions will be obtained numerically by solving for T,B,s, and g in terms of the 
near-horizon parameters (b, q). 



4 Purely Electric and Purely Magnetic solutions 

In this section, we shall briefly review the AdS^ Reissner- Nordstrom (RN) solution for which 
no magnetic field is present, B = 0. For our purposes the RN solution is perhaps better 
referred to as the purely electric solution. We shall then exhibit an exact near-horizon 
solution with B ^ but zero electric charge, whose metric is that of AdS 4 x R 2 . The purely 
magnetic solution is the unique interpolating solution between this near-horizon solution and 
asymptotic AdS6- We shall show, numerically, that the purely magnetic solution exists and 
we shall characterize its asymptotics. 

4.1 The purely electric solution 

In the standard coordinates, the RN solution is given by, 



- 2 + — - — e 2Vl = e 2Vi = r 2 E = — 



In the extremal limit, the location of the horizon r + obeys U(r + ) = U'(r + ) = 0, which 
requires the following relations between q, M and r + , 



q = ^Wrl M = -r 5 + (4.2) 

3 

We see that the extremal entropy density s is proportional to the charge density, s = r+ oc q, 
and does not vanish as T = 0. 



In the coordinates described in section 3.6 the RN solution is given by, 



E = JC^W e 2 *=e 2V > = (Cr + l) 2 (4.3) 
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and U takes the form, 



u = Tc^ + h^^ 2 + &W^ (4 ' 4) 

with 

C = ^r~ Uo = -^ (4 - 5) 

The extremal limit is obtained at q 2 = 10. The near-horizon geometry of the purely electric 
solution is AdS 2 X R 4 . 



4.2 The AdS4 x R 2 near-horizon solution 

Turning on the magnetic field modifies the near-horizon geometry. Simplest is the purely 
magnetic case for which the electric charge vanishes. Here, we can obtain a fairly simple 
picture that will be relevant for interpreting the numerics. Setting Gi — E — 0, a simple 
solution which is regular at the horizon is then given by, 

U = U r 2 e 2Vl =r 2 e 2V * = 1 (4.6) 

where the constants are given by, 




This geometry is recognized as AdS4 x R 2 , and is the D = 6 analog of the AdS3 x R 2 solutions 
studied in [131 EI] • Magnetic brane solutions have been studied further in J22J |23j EU [25], [26] , 
including supersymmetric examples. One lesson is that when these solutions are embedded 
into a higher dimensional supergravity theory they become subject to various potential 
instabilities, and the same could be expected of the AdS 4 x R 2 solution discussed here. 



4.2.1 Perturbing around the AdS 4 x R 2 near-horizon solution 

To construct solutions interpolating between AdS 4 x R 2 and AdS 6 we proceed by including 
small perturbations around AdS 4 x R 2 and using these to seed the flow towards AdSg. Since 
the AdS 4 x R 2 near-horizon solution has E = G\ = 0, the linearized Einstein and gauge field 
equations around this solution decouple from one another. 

We begin by parameterizing the small fluctuations fields for the metric as follows, 

U(r) = U r 2 (l + 2u(r)) 

Vi(r) = ^HU r 2 ) +v l (r) 

V 2 (r) = v 2 (r) (4.8) 
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where the fields u, v±,V2 are treated as first order perturbations. One obtains a closed equa- 
tion for t> 2 from E3, then an equation for V\ from El, and finally an equation for u from 
CON . Putting all together, we find the following system for the Einstein equations, 



„2 // 



r v, 



2 + Arv 2 -9v 2 = 



r 2 v'[ + 2rv[ + r 2 v' 2 ' = 
2ru' + 6u + Arv[ + 6rv' 2 - 6v 2 = (4.9) 

The equations are invariant under dilations in r, as were the original non-linear equations of 



(3.16), and the solutions to the linear equations are linear combinations of powers of r. For 



V2(r) ~ r<J •, we must have 

a 2 + 3a -9 = a = -(-1 + (4.10) 

The positive root for a has been retained here because we are after a regular perturbation 
near r = 0, requiring v% to vanish there. The regular solutions for the other fields are now 
easily found to be, 

u(r) = vi(r) = (cry v 2 (r) = (cr) a (4.11) 

o + 1 

Here, c is an arbitrary integration constant, or scale parameter, which may of course be 
absorbed into the definition of r. Equality of u and V\ guarantees 2+1-dimensional Poincare 
invariances, as expected. 

Since the gauge fields of the AdS4 x R 2 near-horizon solution vanish, the fields E, G\ in the 
gauge field equations may be viewed themselves as parametrizing the first order fluctuations, 
and thus obey the linearized equations, 

(u r 2 E)' + kbd = 

G[+2kbE = (4.12) 



where the constants b and Uq are given by (4.7). Eliminating G%, we find the following 
equation (r 2 E)" — 6k 2 E = for E. This is, of course, again scale invariant, so that its 
solutions are linear combinations of powers of r. The full solution is found to be given by, 

E(r) = r K G x [r) = - ^ r K+1 (4.13) 
where the exponent k must satisfy, 

3 1 

(k + 2)(k + 1) -Qk 2 = K = -- + -Vl + 2Ak 2 (4.14) 
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The positive square root chosen above is dictated by the regularity at the horizon of the 
gauge potential A, which is defined by E = A' together with its vanishing at the horizon. 
Note that we then have k > — 1 for all k, with equality only for k = 0. The electric field 
itself is regular at the horizon provided k > 0, namely when 



k fcjr 




(4.15) 



This condition singles out the value k c for the coupling k (which will turn out to be a critical 
value) as a rough analogue of the special value k = 1 for the AdS 5 case jS]. 



4.3 The purely magnetic brane solution 

We have not succeeded in solving analytically the reduced Einstein equations for the purely 



magnetic brane (the gauge field equations of (3.19) are satisfied automatically for zero charge 



when E = G\ = 0). Numerical analysis is, however, straightforward and shows that the first 
order perturbation of the metric around the AdS4 x R 2 near-horizon solution, computed in 



section 4.2.1 extends into a full fledged solution of the non-linear equations which interpo- 
lates to asymptotically AdSg. Given the above normalizations at the horizon, the AdSg data 
result as follows, 



U(r) 



->■ 1 



n 2Vi(r) 



->• 1 



(4.16) 



and C v is obtained numerically, and given by C v = 1.221905. 



4.4 Entropy density of purely magnetic brane for small T 2 /B 

In the purely magnetic case an expression for the low temperature entropy density may be 
obtained analytically, as we now discuss. This provides a useful check on the numerics. 

The purely magnetic brane solution takes the form, 



ds 2 
F 



dr z 



U 

Bdx 3 A dx 4 



Udt 2 + e 2Vl ((dx 



1\2 



(dx 2 ) 2 ) + e 2V2 



((dx 3 ) 2 + (dx 4 ) 2 ) 



(4.17) 



At zero temperature we have Lorentz invariance in (t, x 1 ^ 2 ), and this fixes e 2Vl = U. The 
zero temperature solution has near-horizon geometry AdS4 x R 2 , so that, as r — > 0, 



U = U r 2 



,2V 2 



-B 



(4.18) 



with Uq = 20/9, as was given in (4.7). 
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At small but finite temperature the near-horizon AdS4 is replaced by the near-horizon 
AdS 4 Schwarzschild solution, given by, 



U = U r 2 



r 



2Vi 



U r 2 



D 2V 2 



-B 



(4.19) 



This is valid provided the horizon lies well within the AdS 4 region, namely for r + C 1. This 
translates into the condition T 2 / B <C 1. The entropy density and temperature are 



Ah 
4 



1 

4 



UoBrlVi 



T 



U'{r + ) 3U r, 



An 



An 



(4.20) 



where V4 is the coordinate volume. Eliminating r + between these expressions, we compute 
the entropy density s as a function of temperature, 



G$S 



10 5 



BT 2 « 1.08 BT 2 



(4.21) 



The T 2 dependence is characteristic of a 2 + 1 dimensional CFT; this CFT is associated to 
the near-horizon AdS 4 region. 



5 Lifshitz solutions for critical k 

The near-horizon Reissner- Nordstrom black brane geometry for k < k c , and the near-horizon 
AdS 4 x R 2 for k > k c are separated by a Lifshitz near-horizon geometry for k = k c . The fields 
for the corresponding near-horizon solution will be computed analytically in this section. 
Linearized fluctuations around this solution smoothly connect onto solutions to the full 
fledged non-linear reduced field equations which interpolate to AdS 6 , as will be shown in 
this section, using the results of Appendix A. 



5.1 Near horizon Lifshitz solution 

In this sub-section we look for scale invariant solutions that can play the role of near-horizon 
geometries in the presence of nonzero electric charge and magnetic field. We start from the 
Ansatz arrived at in section 3J3 for the metric, 

ds 2 = ^~ Udt 2 + e 2V ^({dx 1 ) 2 + (dx 2 ) 2 ) + e 2V *((dx 3 ) 2 + (rfx 4 ) 2 ) (5.1) 
and for the gauge fields, 

F = E(r)dr Adt + b dx 3 A dx 4 

G = Gx(r)dr Adx 1 Adx 2 (5.2) 
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We look for a solution which is invariant under scale transformations of the form, 

11 1 

r -> Ar t -t x 1 ' 2 -> —x 1 ' 2 x 3 ' 4 -> t^x 3 ' (5.3) 

A A Q A^ 

for some real constants a and /?. Scale invariance requires, 

U = U r 2 e 2V 1=r 2a e 2y 2=r 2/3 £ = % Gl = g^a-l (54) 

where we used the freedom to rescale x l to fix the prefactors of e 2Vl ' 2 . The metric is of 
the Lifshitz type [Tj5]. We now substitute these expressions into the reduced field equations 
written in section 13.51 

Assuming bo 7^ 0, equation E2 implies /3 = 0, since all terms in E2 are constant except 
for the final term, which behaves as b"^r~ AI3 . Solving Ml and M2 requires, 

2k 2 b 2 q a 

u = gi = -jT- 5.5 

a kb 



Solving El gives, 



ql = [L-^] k 2 bl (5.6) 



a 



We then find, 

E2-E3 = 2b 2 (3k 2 - 1) (5.7) 
which, assuming as always bo ^ 0, forces, 

k = k c = —= (5.8) 

The remaining equations reduce to the following quadratic equation for a, 

1\ 2 15 

% 

Its solutions may be parameterized as follows, 



a + 2) =M a (5 ' 9) 



, lOz(z-l) , 9 30z 

* = 7IT2F bl = JTTW- (5 ' 10) 



where we have defined, 



z = - (5.11) 

a 

23 



From (5.3) we see that z is the dynamical critical exponent, since under scaling t ~ (x 1,2 ) z . 

We have therefore found a one parameter family of solutions Lif 4 x R 2 parameterized by 
the dynamical exponent z. From (5.10) it is clear that we should restrict to z > 1 in order 
that go and bo be real. All values z > 1 are allowed. For z = 1 we recover the purely magnetic 
AdS 4 x R 2 solution, with (6 , <?o) = (-v/lO/3, 0). As z — >■ oo we go over to the purely electric 

AdS2 x R 4 solution with (b ,q Q ) = (0, vTO). For k k c , only the purely magnetic/electric 
solutions exist. 



5.2 Interpolating solution 

We now look for solutions that interpolate between these k = k c near-horizon solutions and 
an asymptotic AdSe- Here we proceed numerically. We first solve for linearized fluctuations 
around Lif4 x R 2 , identify those that are nonsingular in the near-horizon region, and then 
use these as seeds for the numerical integration out to large r. The details are described in 
Appendix |Aj The upshot is that we find smooth asymptotically AdS 6 solutions for the full 
family of near-horizon solutions described above. 



6 Flows towards low temperature 

We now study the thermodynamics of solutions with nonzero charge and magnetic field. 
In the high temperature limit the behavior is universal, and is governed by the AdS6 
Schwarzschild black brane. The entropy density behaves as s oc T 4 , as appropriate for a 
4+1 dimensional CFT. The low temperature behavior is much more interesting, as it is sen- 
sitive to the charge density, magnetic field, and Chern-Simons coupling k. We have already 
identified k c = 1/a/3 as the critical value of k for which a near-horizon Lifshitz geometry 
arises, and so we can expect a qualitative change in behavior depending on the size of k 
relative to k c , as was summarized in the Introduction. Here we supply more of the details. 



The discussion in sections 3.6 and 3.7 provides the basis for our numerical evaluation of 
the entropy density. As described there, after fixing coordinates the near-horizon data is 
parameterized by specifying the pair (b,q), representing the near- horizon electric and mag- 
netic fields. Numerical integration outwards in r for each choice of (6, q) (and for a given k) 
yields either a smooth asymptotically AdSg solution or a singularity. We discard the singular 
solutions, which requires that we restrict attention to a finite region in the (b, q) plane. For 
each (5, q) within the allowed region we compute the dimensionless quantities (T,s,B,g) 
from the asymptotic data at the AdSg boundary. Our main interest is in understanding the 
behavior of s(T) at fixed B. As we have noted, this depends crucially on the size of k relative 
to k c , and we discuss the three cases in turn. 
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6.1 Low temperature thermodynamics: k > k c 

At vanishing charge density, corresponding to B — > oo, the low temperature entropy density 
was found analytically in (4.21). In terms of dimensionless quantities the result is, 

/ as f ->• at B = oo (6.1) 



3 7T 2 



10 5 

More generally, at finite B the numerics show the low temperature behavior, as T — > 0, 

s = A(k,B)f 2 (6.2) 
Typical examples are shown in figure |4j 



B=0.6 




=0.4 



Figure 4: Entropy density versus temperature at k = 0.7 for selected values of B. 



The T 2 dependence of the entropy density is associated with a near-horizon AdS4 x R 2 



geometry, in particular the one discussed in section 4.2 We note that the electric field 
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vanishes in the near-horizon region, E = 0, which makes it clear that the electric flux 
measured at the AdS@ boundary is sourced by fields living outside the near-horizon region. 
As emphasized in the Introduction, it is the modification of Gauss' law brought about by the 
Chern-Simons term that allows the combined magnetic field and three form field strength 
to act as a charge density. For k > k c , we see that all of the charge density measured at the 
AdSe boundary is generated by this mechanism. 

Another way to exhibit the appearance of AdS4 x R 2 is to plot the flow in (b, q) space 
as we lower the temperature, as shown in figu re |2| Lowering the temperature at fixed B, 
we find that we are driven to the point (6 = J 10/3, q = 0). This point represents the pure 
magnetic brane solution; of course, we never reach precisely this point, since we are keeping 
B fixed at a finite value, while the magnetic brane has B = oo. 

Next we consider decreasing B. Strictly at B = 0, we know that we should recover the 
pure electric RN solution and its nonzero ground state entropy density, and so it is clear that 
A(k, B) must be singular at B = 0. The numerics show the following behavior as B — > 0, 



A(k,B)~ c (k)exp!&> (6.3) 

for some /c-dependent real positive constants c and d; see figure |5j On general grounds we 
expect that this singularity structure should be amenable to analytical understanding, but 
this is lacking at present. 



6.2 Low temperature thermodynamics: k < k c 

For sufficiently small k the effect of the Chern-Simons term is too weak to allow for all the 
charge density to be carried by flux. Some fraction of the flux is instead hidden behind an 
event horizon, and associated with this horizon is a nonzero ground state entropy density. 
Indeed, for k < k c there appears to be such an entropy for any finite B. Representative 
numerical results are shown in figure [2J Note that increasing B suppresses the value of the 
extremal entropy. 

The near- horizon geometry obtained at T = in this regime is AdS2 x R 4 . This can also 
be seen by studying the flow in (6, q) space. As T is lowered at fixed B, we flow towards 
(b = 0, q = a/10), which is the pure electric fixed point. This is illustrated in Figure 



Increasing k at fixed B again reveals k c to represent a critical value. In particular, the 
extremal entropy tends toward zero as k approaches k c from below; this behavior can be 
seen in figure [2] 
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Figure 5: Data illustrating the divergence of A(k, B) as B — >■ 0, in accord with (6.3) 



6.3 Low temperature thermodynamics: k = k c 

At the critical value of k, taking the temperature to zero reveals the appearance of the near- 
horizon Lif4 x R 2 geometry discussed in section |5} Examining the flow in (6, q) space as in 
figure [2] shows that, rather than flowing to either the pure electric or pure magnetic fixed 
point, the IR flow is towards a point on a critical curve. To explain this we first note that 
according to (5.10) the parameters of the T = Lif 4 x R 2 solution obey, 



m + K) 2 -10(6 2 + 3g 2 )=0 



(6.4) 



This same curve governs the T — > behavior of possible values of (b,q), up to the fact 
that there is a finite renormalization relating (b,q) to (&o,9o)- This renormalization comes 
about from the fact that (b, q) are defined at the horizon of a finite temperature solution in 



a coordinate system (defined in section 3.6) that does not agree, in the zero temperature 
limit, with that used to write the zero temperature solution in which (bo, go) appear. The 
coordinate transformation that relates the two solution induces a change in the parameters. 
However, as the temperature is taken to zero, the two coordinate systems differ only in a 
shrinking region near the horizon and, as long as one stays outside this region, the parameters 
(bo,qo) can be read off from the finite temperature solution, and the critical curve ( 6.4[ ) is 
reproduced from the numerics. This procedure is explained in more detail Appendix BJ 



27 




Figure 6: Dynamical critical exponent z = 1/a as a function of B for k = k c . 



The critical curve can be parameterized by the dynamical exponent z, according to (5.10 ). 
The critical curve connects the electric fixed point at b = and the magnetic fixed point at 
q = 0. The full curve is traversed as B ranges from to oo. There is thus a relation between 
the dynamical exponent z and the magnetic field, z = z(B). However, the precise form of 
this relation depends on the full interpolating solution connecting the near-horizon region to 
AdS6- This is because B is defined in terms of asymptotic data, while (6, q) are near-horizon 
data. Numerically, we find that z(B) takes the form shown in figure pi 

The low temperature entropy takes a form dictated by the near-horizon Lifshitz symme- 
try, namely, as T — > we have, 

s = f(B)f 2/z z = z[B) (6.5) 

We have not attempted to characterize the prefactor f(B), but it presumably depends on 
the full interpolating solution. 

These solutions also carry a nonzero value for the three-form charge g. Its behavior is 
shown in figure fj\ The power law dependence for large and small B can be understood as 
follows. First consider the B — > behavior. At B = we have the purely electric RN 
solution, which has G = 0. We can use perturbation theory around the purely electric 
solution to add in a small B. Assuming a finite radius of convergence, we expect that 
the leading behavior of G v is linear in B p , where we now refer to the physical quantities 
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g = -0.68749/5 2 




111(5) 



Figure 7: Plot of the three-form charge g as a function of B at T = and k = k c 



defi ned in (3.27). The p dependence is then fixed by the scale transformations defined 
in (3.28), which give G p ~ p l ^B p . Converting to hatted quantities, this is g ~ l/I? 1 / 2 , 



which agrees with the behavior in figure |7| The same type of argument holds for B — > oo, 
now involving perturbation theory around the purely magnetic solution. The fact that the 
numerics agree with these arguments based on perturbation theory can be taken as evidence 
that perturbation theory is convergent. 



7 Quantum critical behavior 

Our analysis has revealed two interesting examples of quantum critical behavior in this 
system, as manifested by the low temperature thermodynamics. In particular, we focus on 
non-analytic behavior in terms of the control parameters k and B, and the temperature T. 

In the regime k > k c we found nonanalytic behavior at B — 0. This is associated with 
the fact that an infinitesimally small magnetic field is enough to remove the ground state 
entropy density of the pure electric RN black brane. Quantitatively, we saw nonanalytic 
behavior in the coefficient of the T 2 term in the entropy density as B — > 0. As noted already, 
we lack a good analytical understanding of this critical behavior. 



29 



A second example of quantum criticality was found at k — k c , and governed by the 
appearance of a near-horizon Lifshitz geometry. We explored the approach to this critical 
point from several directions, either by decreasing T at fixed B and k, and by letting k 
approach k c at fixed B and (low) T. More generally, there will exist a full scaling form for 
the entropy density in the vicinity of the critical point 

s = F(k-k c ,B,T) (7.1) 

In the analogous AdSs critical point the full scaling function was computed analytically in 
[S], and it would of course be interesting to achieve that for the present case as well, but we 
leave that for the future. 

The two control parameters k and B are on somewhat different footings in the (nominal) 
dual CFT. B is built out of the magnetic field B and charge density p. Changing p cor- 
responds to considering a different state of the CFT, while changing B is accomplished by 
varying an external field. On the other hand, we expect that to change k we need to change 
the field content of the CFT. Chern- Simons terms in the bulk are related to anomalies in 
the dual CFT, and these are determined by the field content. It is therefore unclear whether 
we could ever realize the precise value k c , or whether it is physically sensible to think about 
tuning k. On the other hand, the familiar "fanlike" structure governing finite temperature 
behavior near a quantum critical point shows that it is not necessary to sit precisely at the 
critical point in order to detect its influence. 

The critical D=2+l theories studied here violate parity due to the three- form field 
strength G = G\dr A dx l A dx 2 . In the boundary CFT, this corresponds to an expecta- 
tion value of some two-form operator, (C12) = —(021)- Since we haven't proposed a specific 
CFT dual, we of course cannot say what this operator is. But given such a dual CFT, we 
expect that this two-form operator will play an important role in driving the quantum phase 
transition. 
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A Linearization around near-horizon Lifshitz 

The near-horizon Lifshitz geometry present for k = k c smoothly extends into a solution 
of the full non-linear reduced field equations which is asymptotic to AdS 6 for large r. To 



30 



date, these full non-linear equations have not been solved analytically, but the numerical 
evidence for their existence is overwhelming. The numerical solution is obtained by solving 
the full non-linear reduced equations starting from the near-horizon Lifshitz solution plus 
first order perturbative corrections around Lifshitz. As is familiar from the study of the 
purely magnetic brane in AdSs, the first order fluctuations, subject to appropriate near- 
horizon boundary conditions, will act as seeds for the full non-linear solution with the same 
boundary conditions. 

In this Appendix, we shall compute the first order fluctuations around the Lifshitz solu- 
tion for k = k c . To set up the problem, it is preferable to use the parameter a = z~ l instead 
of the dynamical scaling exponent z. In terms of a, the data of the Lifshitz solution may be 
expressed as follows, 

U(r) = u r 2 E(r) = q 

Vi(r) = a ln(r) G(r) = # r 2a_1 

V 2 (r) = (A.l) 

Here, u ,q ,g , and a are constants related to k and the value b of the field B near the 
horizon with the normalization V^O) = 0. The values of the remaining constants may be 
conveniently parameterized in terms of a, 

V30a yjl0(l-a) 

2a + 1 2a + 1 

20 t—. . 

U ° = (2a + l) 2 g = -\/a(l-a) (A.2) 

Reality of b , q and g clearly requires that a be restricted to the range, 

< a < 1 (A.3) 

This condition restricts 6q to the interval < 6q < 15/4, the upper value being attained when 
a = 1/2. It also restricts q 2 to the interval < q 2 < 10, the upper value corresponding to 
the Reissner-Nordstrom black brane. Recall that the critical curve in terms of near-horizon 
data (bo, qo) is obtained by eliminating a from the above relations, and takes the form, 

3(6^ + qlf- 10(6g + 3g 2 )=0 (A.4) 

A.l Linearized reduced field equations 

The linearization problem is parameterized as follows, 

U(r) = -u r 2 + -ui(r) 
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Vi(r) = ahi(r)+vi(r) 

V 2 (r) = v 2 (r) 

E(r) = g (l + ei(r)) 



G(r) = g y-^l + ^ gi (r)j (A.5) 

Here the fluctuation fields u±, Vi, v 2 , e 1 , and g\ are treated as first order quantities, higher 
orders therein being neglected in the range where all these fields are small. The precise 
normalization of e\ and g\ has been introduced for later convenience. 

The corresponding linearized reduced Einstein field equations are given by, 

El = r 2 v" + r 2 v'z + 2arv[ - 4a(l - a)v! - 2(1 - a)g! 

E2 = 20r 2 v" + 40(2a + l)rv[ + A0arv 2 + a{2a + l) 2 ru[ - 80a(l - a)vi 

+60av 2 + (4a 4 + 12a 3 + 9a 2 + 2a)u x + 10(1 - a)e 1 - 40(1 - a)g x 
E3 = 20r 2 v' 2 ' + 40(1 + a)rv' 2 + 80a(l - a)«i - 180av 2 

+ (4a 4 - 3a 2 - a)u x + 40(1 - a)g x + 10(1 - a)e x 
EA = (2a + l) 2 rV/ + (8a 3 + 24a 2 + 18a + A)ru\ 

+80r(^ + v 2 ) + (8a 4 + 16a 3 + 18a 2 + 10a + 2)m 

+ 160a(l - + 120aw 2 - 60(1 - a)e 1 + 80(1 - a)g Y (A.6) 

The reduced Maxwell-two-form equations are given by, 

Ml = re' 1 + 2r(v' 1 + v 2 ) +Aa(v 1 + v 2 ) + 2ae 1 + 2g 1 
M2 = 20rg[ - a(2a + 1) 2 u[ + A0ar(v[ - v 2 ) 

-(4a 3 + 4a 2 + a)ui + 40a(^i - v 2 ) + 20ae 1 + 20g 1 (A. 7) 

There is also the constraint, 

CON = a(2a + l) 2 ra / 1 +40(l + a)r^+40(2a + l)rw2 

+ (8a 4 + 12a 3 + 6a 2 + a)ui + 80a(l - a)ui - 120aw 2 

+20(1 -a)ei + 40(1 - a)g 1 (A.8) 

whose r-derivative is a linear combination of El, E2, E3, EA, Ml, and M2. 
A. 2 Solution to the linearized equations 

As a result of the scale invariance of the Lifshitz background solution, the linearized equations 
are invariant under an arbitrary rescaling of r. This is manifest also from the observation that 
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the equations above are linear differential equations with constant coefficients with respect 
to the differential operator r d/dr. Thus, the system may be solved by a linear superposition 
of powers of r. Concretely, we set 



(v x {r) \ 
v 2 (r) 
Ui{r) 
e x (r) 

\gi{r)J 



r x V 



V = 



(V X \ 

v 2 



(A.9) 



where Vi,V2,Ui,£i,Gi are independent of r. The allowed values of the exponent A are to 



be determined by substituting (A.9) into equations (A. 6), (A. 7) and (A. 8). The resulting 



equation may be recast in the form AV = 0. The explicit form of A may be deduced from 



(A. 6), (A. 7) and (A. 8), but will not be written out here. Using Maple, one computes, 



det.4. = 400a(2a + 1) 2 Z 2 (/ + 1)(Z + 2a + l) 2 (l 2 + (2a + 1)/ 



6a' 



a 



(A.10) 



For the double zeros of the dete rmin ant, at A = and at A = —1 — 2a, the general solution 
is obtained by replacing r A V in (A.9) by r x V + r A ln(r)V. Here, V and V obey the relations 
AV = and AV + BV = 0, and we shall not need the explicit form of B. 

The modes corresponding to A = and A = — 1 respectively correspond to the dilation 
zero mode of the Lifshitz background, and to its translation mode. The dilation modes are 
non-vanishing at the horizon, therefore modify the horizon data, and are thus excluded from 
contributing. The A = — 1 and A = —1 — 2a modes are both singular at the horizon, and 
cannot contribute either. This leaves the modes for which A obeys, 



I 2 + {2a + 1)/ - 6a 2 



a 







This equation is solved as follows, 



A± = -a - - ± A 



A = ^V28a 2 + 8a + 9 



(A.11) 



(A.12) 



In the interval < a < 1, only the branch l + is regular, since we have Z_ < —2, and l + > 1/2. 
Thus, only the branch l + needs to be retained for matching. The corresponding eigenvector 
is given by the following components, 

V l = -J(A56a 5 + 164a 4 + 78a 3 - 185a 2 - 64a - 44) 

V 2 = J r (192a 5 + 376a 4 + 128a 3 + 106a 2 + 8) + AJ r (12a 4 -8a 3 + 49a 2 + 16a + 12) 
40.7 



a(2a + l) 2 



192a 6 - 560a 5 + 84a 4 - 188a 3 + 75a 2 - 20a + 12 



+A{12a 5 + 76a 4 - 47a 3 - 21a 2 - 16a - 4} 
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4 7 / 

E x = — — 216a 5 - 324a 4 - 70a 3 - 183a 2 - 24a - 20 
2a + IV 

+2A{-36a 4 + 12a 3 - 17a 2 - 4}) 
g 1 = 2J (840a 6 + 244a 5 + 286a 4 - 149a 3 + 12a 2 - 30a + 12) 

+4AJ (l2a 5 + 28a 4 + 21a 3 + 3a 2 + a - 2) (A.13) 



where J is the overall normalization of the eigenvector. 

Maple calculations, using the above regular branch of the perturbative solution as ini- 
tial conditions, clearly demonstrate that the perturbative solution continues into a regular 
asymptotically AdS6 solution. 

B Identification of (60, qo) in low T solutions at k = k c 

Numerical solutions, obtained with initial conditions at the horizon at finite temperature, are 
governed by the data q and b at the horizon, as defined in the coordinate system employed 
in the numerics. These parameters do not coincide exactly with the parameters qo and bo of 
the T = solutions because a finite renormalization effect from finite T boundary conditions 
to the T = scaling solution takes place. This effect is illustrated in the left panel of figure 
[8j At very low T, the functions e 2V2 ^ and E{r) settle to a value which is constant over a 
large number of e-folds, and to high precision. It is the values of e 2V2 ^ and E{r) in this 
regime that determine 60 an d qo in terms of b and q. 

The corresponding assignments for q and b , derived from the numerical runs at finite 
but very low T, may be plotted and compared with the critical curve which is given exactly 



by A. 4 The corresponding plot is given in the right panel of figure pi showing excellent 



agreement throughout the curve. 
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r b 

Figure 8: Left panel: finite renormalization effects occur in transitioning between the finite 
T initial data (b, q) at the horizon, and the initial data (bo, qo) at the horizon of the scaling 
solution. Right panel: the blue line represents the critical curve 3(6q + e^) 2 = lOb^ + 30c 2 ,; 
the data point were collected in the scaling region from the low T numerical solutions. 
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